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ABSTRACT 

Context. Knowledge of the hyperfine structure of molecular lines is useful for estimating reliable column densities 
from observed emission, and essential for the derivation of kinematic information from line profiles. 
Aims. Deuterium bearing molecules are especially useful in this regard, because they are good probes of the 
physical and chemical structure of molecular cloud cores on the verge of star formation. However, the necessary 
spectroscopic data are often missing, especially for molecules which are too unstable for laboratory study. 
Methods. We have observed the ground-state (J = 1-0) rotational transitions of DCO^, HN'^C and DNC with the 
IRAM 30m telescope toward the dark cloud LDN 1512 which has exceptionally narrow lines permitting hyperfine 
splitting to be resolved in part. The measured splittings of 50-300 kHz are used to derive nuclear quadrupole and 
spin-rotation parameters for these species. The measurements are supplemented by high-level quantum-chemical 
calculations using coupled-cluster techniques and large atomic-orbital basis sets. 

Results. We find eQq = -hl51.12 (400) kHz and C, = -1.12 (43) kHz for DCO+, eQq = 272.5 (51) kHz for 
HN"C,andee^(D) = 265.9(83) kHz and e2g(N) = 288.2 (71) kHz for DNC. The numbers for DNC are consistent 
with previous laboratory data, while our constants for DCO^ are somewhat smaller than previous results based 
on astronomical data. For both DCO^ and DNC, our results are more accurate than previous determinations. Our 
results are in good agreement with the corresponding best theoretical estimates, which amount to eQq = 156.0 kHz 
and C, = -0.69 kHz for DCO+, eQq = 279.5 kHz for HN'^C, and eQq(D) = 251.6 kHz and eQq(N) = 309.6 kHz 
for DNC. We also derive updated rotational constants for HN'^'C: B=43545. 6000(47) MHz and £1=93.7(20) kHz. 
Conclusions. The hyperfine splittings of the DCO+, DNC andHN'^^C / = 1 -0 lines range over 0.47-1.28 km s"', 
which is comparable to typical line widths in pre-stellar cores and to systematic gas motions on ~1000 AU scales 
in protostellar cores. We present tabular information to allow inclusion of the hyperfine splitting in astronomical 
data interpretation. The large differences in the '■*N quadrupole parameters of DNC and HN"C have been traced to 
differences in the vibrational corrections caused by significant non-rigidity of these molecules, particularly along 
the bending coordinate. 
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1. Introduction 

Cold interstellar clouds have long been recognized 
as excellent laboratories for determining basic physi- 
cal quantities of molecular structure (for a review see 
iLemaire & CombesI l2007b . In particular, these clouds 
provide access to molecular species that are too unstable 
to permit sufficient terrestrial production for in-depth 
investigation. Some molecules were even detected in 
interstellar space before they were found on Earth. 
A classic example is the X-ogen of the early 1970's, 
which has gained prime importance for studying the in- 
teractions between interstellar g as and magnetic fields 
since its i dentification as HCO+ ( Buhl & Snvdeiil 19701: 
iKramer & Diercksen 1976) . 

Astronomical observations of molecules are not only 
useful to provide accurate rest frequencies of spectral 



Unes dPagani et al.ll200lL l2009h . but also to determine 
hyperfine parameters. The best-known case of hyper- 
fine splitting in molecules with no unpaired electrons is 
the electric quadrupole splitting which occurs for nuclei 
with spin / > 1 such as, for example, D and '"^N. A sec- 
ond type of splitting occurs if the molecule contains nu- 
clei with / > such as H, '^C, and '^N due to magnetic 
spin-rotation coupling and/or nuclear spin-nuclear spin 
coupling. Both effects are usually significantly smaller 
than the splitting due to electron spin-nuclear spin cou- 
pling in molecules with an unpaired electron and an 
/ > nucleus. 

A recent example of astronomically determined hyper- 
fine parameters is the determination of spectroscopic 
parameters of H'^CO^ based on mm-wave observa- 
tions of the dark cloud LDN 1512 (ISchmid-Burgk et al.1 
|2004|) . The exceptionally narrow lines in this cloud 
(AV - 0.16 km s"') allow a frequency resolution and 
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accuracy that usually cannot be attained in the labora- 
tory for short-lived molecules. This category includes 
molecular ions, but also radicals and reactive species 
such as DNC, which has a lifetime of < 1 s in the 
laboratory. Schmid-Burgk et al. showed that even the 
unresolved hyperfine splitting of '^CO plays a role in 
the case of very naiTow lines. They used their spectra 
of LDN 1512 not only to determine the spin-rotation 
constant C/ of the '^C nucleus in H'^CO^, but also to 
show that the magnitude of the splitting and the in- 
tensity ratio of the two resolved features depends on a 
much smaller effect, namely C/ of the H nucleus and the 
nuclea r spin-nuclear spin couplin g between these two 
nuclei jSchmid-Burgk et al.ll2004l) . 



is not particularly accurate, it is in good agreement 
with the very recent, more accurate laboratory value of 
0.2645 (46) MHz (g echtel et al. 2006). More recently, 
ICaselh et al.l (Il995h determined hyperfine parameters 
for the N2H^ J = 1-0 line from astronomical data, 
which iGerin et all ( 200 Ih e xtended to higher-7 transi- 
tions, and lDore et al.l (12004) to the N2D+ isotopologue. 
This paper describes a new determination of the hyper- 
fine structure of the 7 = 1 - lines of DCO+, DNC and 
HN'^'C, based on observations of LDN 1512. The mea- 
surements are supplemented by high-level quantum- 
chemical calculations of the corresponding hyperfine 
parameters using state-of-the-art coupled-cluster tech- 
niques together with large atomic-orbital basis sets. 



In astrophysics, the prime use of hyperfine splitting 
is the possibility of measuring the optical depths of 
molecular lines. This fundamental quantity allows es- 
timation of molecular column densities without as- 
sumptions about the beam filling factor. The classic 
example is the split ting of the NH3 inversion lines 
(iHo & TownesI [T983I) which are used widely to mea- 
sure the kinetic temperatu r es of d ense interstellar clouds 
([Walmsle v & UngerechtsI Il983l) . In addition, knowl- 
edge of hyperfine structure is essential for deriving 
the kinematic structure of clouds from observations of 
molecular line profiles. In particular, the central re- 
gions of pre-stellar cores are cuiTently of great inter- 
est, as the places where the transiti on from spheri- 
cal in fall to disk-like rotation occurs ( Bergin & Tafallal 
l2007h . However, in these objects, many 'standard' kine- 
matic probes are unavailable du e to freeze-out onto 
dust grains (i Bergin et al.1 12002|) . Deuterium bearing 
molecules are abundant even under these conditions, but 
their lines always exhibit hyperfine splitting due to the 
nonzero spin (7 = 1) o f the D nucleus. For example, 
IVan der Tak et al.l (l2005h used observations of H2D^ to 
study the kinematics at the center of the pre-stellar core 
LDN 1544, where most other molecules are frozen onto 
dust grains. 

In the laboratory, the Doppler-limited line widths in 
the 3 mm region are of the order of 200 kHz, mak- 
ing it possible to resolve moderate to large hyper- 
fine splitting. However, very small quadrupole splitting, 
such as the '"^N splitting in HNC or essentially all of 
the deuterium quadrupole splitting cannot be resolved 
in this frequency region. Sub-Doppler resolution tech- 
niques, such as Lamb-dip spectroscopy or molecular 
beam millimeter-wave Fourier transform spectroscopy, 
are sometimes available, but the former is usually not 
feasible for short-lived species, and the latter type of 
measurements is available only in very few laboratories. 

While the relatively large nitrogen quadrupole split- 
ting of DCN has been resolved in the laboratory 
(iBriinken et al.l 12004) and in space (iTurner, 2001). the 
hyperfine structure of the astrophysically important 
DCO+, HN'^C and DNC species had not yet been re- 
solved at the time of our observations (2004-2005). 
High-resolution spectroscopy of these species is not 
only astronomically useful to determine the optical 
depths of the lines, but also to verify the results of 
quantum-chemical calculati ons of spectroscopic con- 
stants. For example, Frerkin Fet al.l ([1979) measured a 
nitrogen quadrupole moment of 0.28 (3) MHz for HNC 
in the dark cloud LDN 134N. Even though this value 



2. Observations 

Observations of the 7 = 1-0 lines of DCO+, DNC 
and HN'^C near 72039, 76306 and 87091 MHz were 
performed on 2004 August 5-6, and 2005 March 12 
and July 29. The 30-m telescope of the Institut de Radio 
Astronomic Millimetrique (IRAM£] was used, with the 
facility receivers A 100 and B 100 as front end. The tun- 
ing range of these receivers had just been extended from 
80 GHz down to 70 GHz. The beam size is 33" in this 
frequency range, and the main beam efficiency is ai80%. 
We used the Versatile Spectral Assembly (VESR\) cor- 
relator as backend to achieve a spectral resolution of 
3.3 kHz, or 0.013 km s"'. Sideband gain ratios, mea- 
sured with the Martin-Puplett interferometer, were in 
the range 0.91-1.13, depending on backend module. 
System temperatures were s;300 K at 72 GHz, x;200 K 
at 76 GHz and x;150 K at 87 GHz. Integration times 
were 78 minutes (on-i-off) at 72 GHz, 215 minutes at 
76 GHz and 197 minutes at 87 GHz, giving rms noise 
levels of rMB = 158 mK at 72 GHz, 47 mK at 76 GHz 
and 41 mK at 87 GHz. Pointing was checked every hour 
on the nearby planet Venus. 

The position observed is or = 05:00:54.40, 6 - 
H-32:39:37.0 (B1950). The cloud velocity at this posi- 
tion, as determined from observations of the HC3N 3-2 
line near 27294 MHz with the Effelsberg 100-m tele- 
scopeQ, is Vlsr= 7.069+0.001 km s"'. The beam size 
of the 100-m telescope at the HC3N frequency is almost 
equal to that of the 30-m telescope at the DCO+, DCN 
and DNC frequencies, so the effect of the known ve- 
locity gradients within the LDN 1512 cloud should be 
very small. An upper limit is obtained by convolving 
the HC3N data to a 50" beam, which gives a velocity of 
Vlsr = 7.063 + 0.001 km s"' . Therefore the effect of the 
velocity gradient in LDN 15 12 is less than 0.006 km s ' . 
For further discussion of these ve locity gradients and 
for acc urate HC3N frequencies, see lSchmid-Burgk et al.l 
(12004 . 

The data reduction was performed using the Continuum 
and Line Analysis Single-dish Software (CLASS) pack- 



' IRAM is an international institute for research in mil- 
limeter astronomy, co-funded by the Centre National de la 
Recherche Scientifique (France), the Max Planck Gesellschaft 
(Germany) and the Instituto Geografico Nacional (Spain). 

- The Effelsberg telescope is operated by the Max-Planck- 
Institut fiir Radioastronomie 
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ag43 and followed standard procedures, except that the 
combination of data taken on different epochs posed un- 
usual challenges. The telescope system corrects the fre- 
quency scale of the spectra for heliocentric motion of 
the Earth, but the correction is only exact for the center 
of the receiver bandpass. In our case, the HN'^^C line 
was observed at an offset of 160 MHz to allow simul- 
taneous observation of other lines. The telescope sys- 
tem adds this offset as an absolute value, but the or- 
bital motion of the Earth actually dilates or shrinks the 
frequency scale by a factor of ~10"'*, the ratio of the 
telluric orbital velocity and the velocity of light. The 
HN'^C spectra were taken 121 days apart, and the fre- 
quency shift is 16 kHz x cos 5, where 6 is the eclip- 
tic latitude of the source (10°). Indeed, the two HN'^C 
spectra agree very well after a shift of 21 kHz, and sim- 
ilarly, the DNC spectra are consistent with the expected 
shift of 4.4kHz for an epoch separation of 340 days. To 
obtain the correct absolute frequency scale, the spectra 
of all epochs were shifted to a hypothetical observing 
date of June 8, which is the day in the year that the Sun 
and LDN 1512 have the same hour angle and transit at 
the same time. We assume the frequency accuracy to be 
3 kHz. This is about twice the LSR velocity gradient of 
the source and leaves some room for possible errors in 
the correction of the hehocentric motion. 



fields. The latter are o btained using analytic second- 
derivative techniques (Gauss & Stanton 1997j as de- 
scribed in Stanton & Gauss (2000). Spahn et al] (|2008|) 
provide recent examples of the very good agreement 
between experiment and calculations for DC3N and 
HC3N. 

The nuclear spin-rotation tensor is a second-order re- 
sponse property and can be computed via the associated 
second derivative of the energy with respect to nuclear 
spin and the rotational angul ar mo mentum as perturba- 
tions. As described in lGauss et al.l (Il996l) . we calculate 
the spin-rotation tensor using perturbation-dependent 
basis functions (so-called rotational London atomic or- 
bitals) in order to ensure fast basis-set convergence. 
No vibrational averaging has been performed for the 
spin-rotation tensors, as the corresponding theoretical 
expressions exhibit numerical problems in the case of 
polyatomic linear molecules. 

To ensure convergence in the electron-correlation treat- 
ment, a few calculations have been performed at 
levels be yond CCSD( T), i.e. at the CCSDT level 
jNoga & B artlett 1987; ,Scuseria & Schaefed 19881) with 
a full treatment of triple excitations and at th e CCSDTQ 
level (iKallav & Suriinl2001l:lK'allav et al.l2004) with an 
additional consideration of quadruple excitations. 



1998 



2007) 



3. Quantum-chemical calculations 

High-level quantum-chemical calculations using 
coupl e d-cluster (CC) techniqu es (Cizek 1966a b; Gauss 
ICrawford & Schaefed 12000; .Bartlett & Musiat 
have been performed for the hyperfine pa- 
rameters of DCO+, HN'^C, and DNC. Calculations 
have been carried out for the mos t part at the CC 
single s and doubles (CCSD) level dPurvis & Bartlett! 
Il982l) augmented by a pert urbative treatm ent of triple 
excitations (CCSD(T); Ragh avachari et alJll 989) which 
has proven in many cases to provide a reliable account 
of electron-correlation effects on energies and prop- 
erties. The required one-particle basis sets for these 
calculations have been taken from Dunning's hierarchy 
of correlation-consistent basis s ets with cc-pV XZ 
denoting the standard valence sets dDunnin j|l98 9), cc- 
pCVXZ those with additional core-polarizing functions 
(IWoon & DunningI 1 1 995h . and au g-cc-p(C)VXZ thos e 
with additional diffuse functions dKendall et a n il992h . 
Here, X represents the cardinal number of the basis sets 
and values of 3 (= T), 4 {- Q), and 5 have been chosen 
in the present work. 

The theoretical determination of quadrupole coupling 
constants is based on the evaluation of the electric- 
field gradients (efgs) at the corresponding nuclei. The 
nuclear quadrupole moments required to convert the 
efgs to the quadrupole couplings are taken from the 
literature dPvvkkol 120081) . and the following values 
have been here adopted: 2.860 (15) mbarn for D and 
20.44 (3) mbarn for ^'^N. The reported calculations also 
take into account zero-point vibrational effects on the 
quadrupole couplings. Those are treated in a perturba- 
tive manner as described in lAuer et alj d2003) and ne- 
cessitate the evaluation of quadratic and cubic force 



All calculations have been carried out using the CFour 
quantum-chemical program packag^ Only the CCSDT 
and CCSDTQ calculations have been carried out with 
the MRCC packag^ which has been interfaced to 
CFouR. 

Tables [1] and |2] summarize our computational results for 
the hyperfine parameters of DCO+, HN'^^C, and DCN 
and in particular document the convergence of the cal- 
culated values with respect to basis set and electron- 
correlation treatment. 



4. Results 

4. 1 . Hyperfine splitting of DCO^ 

Figure[T]shows the DCO^ spectrum toward LDN 1512. 
The spin / = 1 of the D nucleus splits the transi- 
tion into three components, labeled by the total angu- 
lar momentum F = J + I. The three expected hyper- 



http://www.iram.fr/IRAMFR/GILDAS 



Coupled Cluster techniques for Computational 
Chemistry, a quantum-chemical program package by J. 
F. Stanton, J. Gauss, M. E. Harding, and P. G. Szalay with 
contributions from A. A. Auer, R. J. Bartlett, U. Benedikt, 
D. E. Bemholdt, C. Berger, Y.J. Bomble, O. Christiansen, 
M. Heckert, O. Heun, C. Huber, D. Jonsson, J. Juselius, K. 
Klein, W. J. Lauderdale, D. Matthews, T. Metzroth, D. R 
O'Neill, D. R. Price, E. Prochnow, K. Ruud, F. Schiffmann, S. 
Stopkowicz, M. E. Vamer, J. Vazquez, J. D. Watts, F. Wang 
and the integral packages MOLECULE (J. Almlof and P. R. 
Taylor), PROPS (P R. Taylor), ABACUS (T. Helgaker, H. J. 
Aa. Jensen, P. J0rgensen, and J. Olsen), and ECP routines by 
A. V. Mitin and C. van Wiillen. For the current version, see 
http : / /www . cf our . de 

' MRCC, a string-based quantum chemical program suite 
written by M. Kallay. See http : //www . mrcc . hu 
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fine components are clearly detected and resolved. We 
derive the frequencies of these components by fitting 
three independent Gaussians to the observed spectrum. 
The fit uses a cloud velocity of Vlsr = 7.069 km s"' 
(see § |2]l and a line width of 35 kHz which is the av- 
erage width of the three components. The spacings of 
the components, independent of Vlsr, are -0.198 (6) 
and +0.269 (6) km s"' which is +47.58 (144) and 
-64.64 (144) kHz (Table|3]l. These splittings are Hnear 
combinations of the quadrupole coupling constant eQq 
and the spin-rotation cons tant C/ of the D nucleus (e.g., 
iTownes & Schawlowlll955l) . Solving the two equations 
with two unknowns gives the eQq and C/ values given 
in Table S 

The present value of C/ in Table B] is srn aller in mag- 
nitude than that bv lCaselli & Dorj72005l) while eQq is 
larger. However, agreement exis t s with in the larger un- 
certainties from ICaselU & Dorg (l2005h . These authors 
also observed the LDN 1512 cloud with the IRAM 30m 
telescope, but at a different position, where the lines 
are weaker. Therefore, the signal-to-noise ratio of their 
spectra is lower, even though the integration time has 
been longer. 

The quantum-chemical calculations for DCO^ (see 
Table [U yield as best values for eQq 156.0 kHz 
(CCSD(T)/cc-pCV5Z calculations plus vibrational cor- 
rections); the small changes from the quadruple-zeta 
(QZ) to the quintuple-zeta (5Z) basis set suggests a 
minute increase of probably less than 1 kHz upon 
extrapolation to the basis-set limit. The agreement is 
good, within twice the experimental uncertainties, for 
the fit in which C/ has been varied. The agreement is 
slightly worse for the fit in which C/ was kept fixed to 
the quantum-chemically calculated value and slightly 
worse still in compar ison to the previous values from 
ICaseni&Dorg(l2005h . see TableEl 

The best theoretical value for C/ is -0.69 kHz. As dis- 
cussed in § [3] no vibrational corrections could be com- 
puted for this value. However, since the vibrational cor- 
rection to eQq is fairly small, it is likely that the one 
on C/ will be small also, possibly of order of 0.1 kHz. 
The convergence patterns seen in Table |2] suggest that 
remaining errors in the calculations can be considered 
rather small so that we conclude, based on these theo- 
retical esti mates as well a s the re sults of our measure- 
ments that lCaselli & Dor^ ( l2005h overestimated C/ due 
to the limited signal-to-noise in their spectrum. Both 
fitted values agree within the fairly large uncertainties 
with the theoretical value. 

The absolute positions of the strongest hyperfine com- 
ponent from the present astronomical observations 
agree to about 3 kH z with the previous observation 
(ICaselh & Dorell2005h . 



4.2. Anomalous excitation of DCO"^ 

The continuous line in Figure[T]is a model fit to the data 
that uses the splitting constants that we just derived. The 
model was computed with the hyperfine structure (HFS) 
routine within CLASS. This routine assumes that the 



hyperfine components of the transition have the same 
width and the same excitation temperature. While the 
data support the first assumption, the second appar- 
ently does not hold. The highest-frequency F = I - I 
component has a higher intensity than the model pre- 
dicts. Such 'hyperfine anomalies' are well-known for 
the /= 1 - line of HCN, where the F = - 1 
component is unexpectedly strong toward many dark 
clouds (e.g.. lWalmslev et aLlll982l) . The implied differ- 
ence in excitation temperatures between the hyperfine 
components can be understood through detailed cal- 
culation s of hyperfine selective collisional cross sec- 
tions (Monteiro & Stutzkilll986l) . Such calculations do 
not exist for DCO^ at present, although they have been 
reported for NaH^ (iDaniel et al.l 12005 "*). where a simi- 
lar ex citation anomaly has been observed ( Case lli et al.l 
119951) . Hyperfine selective collision data for the electron 
impact excit ation of HCN, HN C, DCN and DNC are 
presented bv lFaure et al.l (l2007h . but do not apply here, 
since the electron abundance in dark clouds is low. 



4.3. Hyperfine splitting of DNC 

Figure |2] shows the spectrum of DNC toward LDN 
1512. The profile consists of one weak, isolated compo- 
nent at low frequency, and two stronger, broader peaks 
which appear to be blends of several components. The 
spins / = 1 of the D and '"^N nuclei lead theoretically to 
a splitting into nine components at seven different fre- 
quencies, but not all these components are resolved in 
the current data. Since the hyperfine splitting of both 
nuclei are of similar magnitude, we have used the sym- 
metric cou pling scheme /tot = fy + Id and F - J + I^ot- 
In contrast. 'B echtel etTa D (12006^ use the sequential cou- 
pling scheme, which only affects the labeling of the 
levels, but not the shape of the hfs pattern. The iso- 
lated component is well fitted with a Gaussian profile of 
width 0.150 (1) km s"'. We have fitted the profile with 
an 'empirical' model consisting of several Gaussians 
with independent positions and intensities, but with the 
width fixed to 0. 15 km s"' . The fit has significant resid- 
uals if four or five Gaussians are used, but the sum of 
six Gaussians (plus one for the isolated feature) gives a 
satisfactory fit to the observed profile (Fig.|2^). 

The results of the empirical model are in good agree- 
ment with our quantum-chemical calculations, which 
predict seven hyperfine components for this transition, 
due to the nuclear spin of both D and N. In the case 
of DNC, the quadrupole moments eQq of the D and N 
nuclei are similar in magnitude, unlike in DCN where 
eQq(N) » eQq{D). 

The positions of the components of the empirical model 
are in good agreement with the predictions based on 
quantum-chemical calculations, except for the highest- 
frequency component. Suspecting that this mismatch 
is caused by the small separation of the 2,2^2,2 and 
the (coinciding) 0,1^0,0 and 0,1— >2,2 components, we 
have reduced the empirical model to five Gaussians, and 
freed up the width of the fourth Gaussian, which rep- 
resents these components. This 'semi-empirical' model 
gives an excellent match to the data, as shown by the 
solid line in Fig. |2t). The width of the fourth Gaussian 
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is 0.202 (9) km s"' implies a separation of the 2,2^2,2 
component from the 0, 1 — >0,0 and 0, 1 — >2,2 components 
of 0.05 km s"' . This result is in good agreement with the 
theoretical results: the calculated splitting is 8.7 kHz, 
but the apparent separation is smaller because of the 
small intensity of the 0,1— >0,0 and 0,1— >2, 2 compo- 
nents; see Table [3] 

The hyperfine parameters determined for DNC are com- 
pared in Table |5] with those previously determined in 
the laboratory by Bechtel et al. (2006) and with the 
best ab initio values. We note a good agreement be- 
tween the two experimental sets of data, but emphasize 
that our values have reduced error bars. The compar- 
ison with the theoretical best estimates is also favor- 
able and lends support to the reliability of the exper- 
imentally determined hyperfine parameters. The theo- 
retical quadrupole coupling terms determined with the 
largest basis set and incorporating vibrational correc- 
tions are 257.6 and 309.6 kHz for the D and '"^N nu- 
cleus, respectively, while the calculated nitrogen spin- 
rotation constant amounts to 5.54 kHz. The theoreti- 
cal valu es are appar ently slightly closer to the previous 
values jBechtelet a l. 2006). While a smooth and fast 
convergence is observed in the calculation of the spin- 
rotation and deuterium quadrupole coupling constants, 
the determination of the nitrogen quadrupole coupling 
constant turned out to be challenging. Going from the 
quadruple-zeta to the quintuple-zeta set causes a change 
of 16 kHz and vibrational corrections amount to more 
than 40 kHz, while electron-correlation effects beyond 
CCSD(T) are found to be less important. However, the 
noted problems in the accurate determination of the ni- 
trogen quadrupole coupling constant can be traced back 
to the smallness of the actual value which is roughly 
one order of magnitude smaller than the correspond - 
ing quadrupole coupling in DCN (iBriinken et al.ll2004h . 
Moreover, the changes in the calculated eQq values de- 
crease fairly rapidly by a factor of about 6 with increas- 
ing size of the basis set. Taking into account the de- 
crease in the magnitude of the vibrational corrections 
one can estimate a value of about 306.5 + 1 kHz as 
the CCSD(T) value at the inifinitely large basis set, in 
agreement within the uncertainty of the laboratory value 
and within 2.5 times the uncertainties of the present as- 
tronomical observations. 

The absolute positions of the isolated hyperfine com- 
ponents from the present astronomical observations 
agre e to 1 kHz with th e previous laboratory frequen- 
cies (iBechtel et al.l2006h . The laboratory frequency was 
deemed to be uncertain to 5 kHz while our astronomical 
frequency is uncertain to probably no more than 3 kHz. 



4.4. Hyperfine splitting of HN^^C 

Figure |3] shows the spectrum of HN'-^C 7=1-0 
toward LDN 1512. One weak, isolated hyperfine fea- 
ture at low frequency (LF for short) as well as two 
somewhat resolved ones can be recognized easily as 
would be expected for a molecule with a '^N nu- 
cleus. Analyses of the astronomical spectrum as well 
as simulations of the spectrum based on t he experimen- 
tal quadrupole moment {eQq) for HNC (iBechtel et al.l 



l2006l) . the quantum-chemically calculated nuclear spin- 
rotation parameters (C/) and the spin-spin coupling 
terms (Sxy) derived from the interatomic distances and 
the nuclear magnetic moments reveal that the situation 
is more complicated. First, contributions from the '^C 
nucleus considerably broaden the central hyperfine fea- 
ture (CF for short) such that a shoulder to the lower 
frequency side should be discernable. Second, the hy- 
perfine pattern should get wider, and the relative in- 
tensity of the highest frequency feature (HF for short) 
should increase somewhat with respect to the central 
one. Splitting caused by the H nucleus increases the 
total number of hyperfine components further, but the 
total number of observable features does not increase. 
However, it is noteworthy that the H hyperfine parame- 
ters counteract somewhat the effects caused by the '■'C 
nucleus mentioned under "second". The net effect is one 
isolated, weak hyperfine feature at low frequencies, one 
strong and very broad CF with a low-frequency shoul- 
der, and a HF that should be stronger than the isolated 
component, but weaker than the central one. 

Analyses were started assuming three different hyper- 
fine features with three different widths. The LF-HF 
splitting seemed reasonable, but the HF was too strong 
compared with the CF, and the position of the CF was 
too low in frequency. Fitting two features for the CF re- 
quired their width to be equal to that of the LF which 
was determined to be 43 kHz as expected from the sim- 
ulations. A fit of four components with widths of 43 kHz 
was still fairly poor. Simulations suggested the width of 
the HF to be about 15 % larger than that of the remain- 
ing components. The analysis of the astronomical spec- 
trum improved somewhat, but was still not satisfactory. 
Moreover, a fit of eQq to the splittings with all other hy- 
perfine parameters kept fixed yielded a value of 272.5 + 
5.1 kHz, which is larger than the one d etermined for 
HNC, 264.5 + 4.6 kHz, in the laboratory teechtel et aTl 
2006), but in agreement within the quoted uncertainties. 
The rms of the fit is 5.1 kHz. Releasing the width of 
the HF gave the best analysis of the astronomical spec- 
trum, but yielded a width of 7 1 .5 +2.3 kHz, much bigger 
than expected from the simulations. The spectroscopic 
fit yields eQq - 258.7 + 5.1 kHz, which is now smaller 
than the HNC value, but again in agreement within un- 
certainties. The rms of the fit is, with 4.0 kHz, slightly 
better than the fit mentioned before (Table |3]l. 

The contributions to the hyperfine splittings decrease 
from ''^N over '^C to H. Therefore, we apply a sequen- 
tial coupling scheme: Fx^J + /("*N); F2 = J + iC^C); 
F - J + /(H). Only eQq could be determined for the 
i^N nucleus. Trial fits with C("'N) or C(i^C) or both 
released yielded values for these parameters that were 
deemed to be unreliable as they could differ by as much 
as about 5 kHz from the predicted values. 

The final values for the hyperfine parameters for HN'^C 
are given in Table |6] The resulting nitrogen quadrupole 
coupling of 272.5 (51) kHz is in good agreement with 
the best theoretical value of 279.5 kHz, in particular 
when noting the slow basis-set convergence in the corre- 
sponding calculations. The vibrational contributions are 
much larger still than the akeady large contributions for 
DNC; they amount to about 46 kHz for HN'^C. These 
contributions decrease again slightly in magnitude go- 
ing from the quadruple-zeta to the quintuple-zeta basis 
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set. Thus, the CCSD(T) eQqO*N) value at the basis-set 
Umit is probably around 3 kHz smaller than the one cal- 
culated at the cc-pCV5Z basis set with vibrational cor- 
rections and now agreeing within the uncertainty with 
the value from astronomical observations. 



4.5. Rotational and distortion parameters for H/V'^C 

Accurate rotational as well a s distortion parameters are 
availa ble for DNC (Brunke n et al l 120061 i Bechtel et al 



20061) and DCO+ (CaselH & Doreii2005riLattanzi et al 



2007h and the present measurements will only mod- 
ify these parameters slightl y. In contrast, only sparse 



data are available for HN'^C. IPearson et alj ( 119760 mea- 
sured the y= 1 <— and 2<— 1 rotational tr ansitions of sev- 
eral HNC isotopologs. More recently, iMaki & Mellaul 
([2001) obtained rovibrational transitions of the V2 bend- 
ing mode of HN'^C together with extensive data for 
HNC. Replacing the J-l(—Q laboratory rest frequency 
with our astronomical rest frequency (Table [3]) reduces 
the uncertainties of B and D by factors of 5 and 1.7, 
respectively . Fixin g // to a value of 157 mHz take n 
from HNC dThorwi rth et al. 2000; Bechtel etaDHoO^, 
we obtain B = 43545.6000 (47) MHz and D = 
93.7 (20)kHz; the numbers in parentheses are one stan- 
dard deviation in units of the least significant figures. 
The present values of B and D agree with the previ- 
ous o nes within their larger uncertainties dPearson et al.l 
11976^ The resulting predictions are accurate up to 
ss500 GHz but should be viewed with some caution at 
higher frequencies. 



5. Discussion and conclusions 



5.1. Astronomical implications 

At spectral resolutions of ~0.1km s as commonly 
achieved by astronomical instrumentation in the 3 mm 
wavelength range, the J - 1-0 line of DCO^ splits into 
three hyperfine components, the HN'^'C 7=1-0 line 
into four, and the DNC J = 1-0 line into six. Table [T] 
gives the velocity offsets and the relative strengths of 
these components. The strengths are theoretical values 
which only apply if the lines are optically thin and the 
excitation is thermalized; as Section l4~2l discusses. these 
conditions are often not fulfilled, so that detailed radia- 
tive transfer modeling is needed to extract physical in- 
formation out of the observed hyperfine intensities. 

The hyperfine splitting extends over 0.47 km s ' for 
DCO+, over 1 .28 km s^^ for DNC, and over 0.73 km s ' 
for HN'^C. These intervals are comparable to or larger 
than the total (thermal -i- turbulent) line widths in many 
astrophysical objects, in partic ular pre- and protostellar 
cores in star-forming regions dBergin & Ta falla 2007). 
Taking these splittings into account is therefore crucial 
to derive accurate column densities of DCO^, DNC and 
HN'^C from their ground-state rotational lines. 

The magnitude of the hyperfine splittings is also com- 
parable to the velocities of infalling, outflowing and ro- 
tational motions around protostars on ~1000 AU scales 



(e.g.. lHogerheiidell2001h . Taking the splitting into ac- 
count is therefore essential for a correct determina- 
tion of the velocity field of the gas from these lines. 
The quadrupole splittings of linear molecules decrease 
rapidly with increasing J, at least for the strong tran- 
sitions, so that going to higher-/ lines would seem a 
way to avoid this complication. However, at the typical 
(low) temperatures and densities of pre- and protostellar 
cores, the higher-7 Hues of DCO+, DNC and HN'^C are 
usually less excited and thus weaker than the ground- 
state lines. 



5.2. Molecular physics implications 

The most remarkable aspect of the present results from 
the molecular physics point of view are the large vi- 
brational coiTections to the '"*N eQq values of HN'-'C 
and DNC and the corresponding large differences in the 
ground state eQq values. Isotopic differences in experi- 
mentally determined quadrupole coupling constants are 
usually not or barely significant. Significant effects have 
been observed e.g. by Gatehouse et al. (1998) where 
'"B/"B isotopic shifts of around 15 kHz were deter- 
mined for quadrupole coupling constants of both ^^Cl 
and ^^Cl of CIBO, though these shifts were only very 
small fractions (< 10"^) of the constants themselves. 
The authors attributed these effects to the non-rigidity of 
the molecule in particu lar along the bendi ng coordinate. 
Experimental work (iBechtel et al.l |200"8|) and calcula- 
tions (Wong 2008) lend support to this view and sug- 
gest that low-order vibrational corrections to the eQq 
values are probably not sufficient for very accurate pre- 
dictions of quadrupole coupling constants for specific 
vibrational states. 

While the bending mode of HNC is known to be fairly 
non-rigid, those of the isoelectronic HCN, HCO^, and 
NaH^ molecules are much more rigid. However, the 
bending mode of the HOC^ molecule is even less rigid 
than that of HNC. Therefore, it would be interesting 
to investigate the '^O quadrupole coupling constants 
of Hi^OC+ and Di^OC+. The HOC+ molecule is less 
abundant than its isomer HCO^ by factors of ~100 in 
diffuse clouds and ~1000 in dense clouds (iLiszt et al.l 
,2004), but HOC^ is strongly enhanc ed in regions with 
strong radiation fields an d/or shocks (iRizzo et alj|2003l; 
iFuente et al.ll2005l 12008 ). so these may be good places 
to search for its rare isotopologues. In addition, this 
molecule is not so easy to produce in the laboratory. 
Thus, the easiest way to obtain information on the 
hyperfine structure of this ion are detailed quantum- 
chemical calculations. 



New or updated predictions of the rotational spectra for 
DCO+, DNC, and HN'^C will be available in the cata- 
log section of the Cologne Database for Molecular spec- 
troscop>0 dMiiller et alJl200lLl2005h . 
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Fig. 1. Spectrum of the 7 = 1-0 transition of DCO^ 
toward the dark cloud LDN 1512, obtained with the 
IRAM 30m telescope. Vertical scale is antenna tem- 
perature (Tmb) in K. Overplotted is a model fit which 
assumes that the three hyperfine components have the 
same width and excitation temperature. 
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Fig. 3. Spectrum of the 7 = 1-0 transition of HN'^C 
toward the dark cloud LDN 1512, obtained with the 
IRAM 30m telescope. Vertical scale is antenna temper- 
ature (Taib) in K. 
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Table 1. Computed quadrupole-coupling constants eQq (kHz) for DCO"", HN"C, and DNC.° 



Computational level 


DCO* 


HN"C 


DNC 




D 


N 


D 


N 


CCSD/cc-pCVDZ 






J\JO.\J 




CCSD(T)/cc-pCVDZ 






309.3 


1021.4 


CCSDT/cc-pCVDZ 






309.3 


1029.1 


CCSDTQ/cc-pCVDZ 






309.2 


1033.2 


CCSD(T)/cc-pCVTZ 






290.8 


431.9 


CCSD(T)/cc-pCVQZ 


159.2 


341.6 


279.8 


341.6 


CCSD(T)/cc-pCVQZ + vib 


153.6 


294.8 


259.1 


325.2 


CCSD(T)/cc-pCV5Z 


162.2 


325.5 


278.3 


325.5 


CCSD(T)/cc-pCV5Z + vib 


156.0 


279.5 


257.6 


309.6 



" All calculations have been performed at geometries obtained at the same computational level. 
Table 2. Computed nuclear spin-rotation constants d (kHz) for DCO+, HN^^C, and DNC." 



Computational level 


DCO+ 
D 




HN'^'C 




DNC 


H 


N 


C 


D 


N 


CCSD(T)/cc-pCVTZ 


-0.701 


-4.551 


22.954 


6.194 


-0.611 


5.417 


CCSD(T)/cc-pCVQZ 


-0.695 


-4.505 


23.345 


6.292 


-0.605 


5.502 


CCSD(T)/cc-pCV5Z 


-0.693 


-4.486 


23.471 


6.330 


-0.602 


5.536 


CCSD(T)/aug-cc-pVTZ 


-0.700 


-4.575 


22.878 


6.200 


-0.614 


5.423 


CCSD(T)/aug-cc-pVQZ 


-0.695 


-4.510 


23.318 


6.295 


-0.605 


5.505 


CCSD(T)/aug-cc-pV5Z 


-0.693 


-4.486 


23.456 


6.331 


-0.602 


5.537 



° The calculations have been performed at geometries obtained at the CCSD(T)/cc-pCVQZ level of theory (DCO* : r(CO) = 
1.10639 A, r(HC) = 1.09236 A; HN"C, DNC: r(NC) = 1.16927 A, r(HN) = 0.99526 A). 

Table 3. Measured hyperflne splitting A" (kHz) and residual o-c of DCO"^, DNC, and HN'^C relative to a selected 

hyperflne component. If no numbers are given for the value of a measured splitting and its residual, the hyperflne 
component overlaps with the previous component. The values of the splitting and the residual refer to the respective 
intensity weighted averages. 



Transition 


Splitting 


Residual* 


DCO+: F = 2- 


1 at 72039.306 (3) MHz 


1 - 1 


+47.58 (144) 


1.20 


0- 1 


-64.64(144) 


0.80 


DNC: l,F = 2, 


1 - 2, 2/2, 1 - 0,0 at 76305.512 (3) MHz 


1, 1 - 1, 1 


115.30(100) 


-0.18 


2, 3 - 2, 2 


175.44 (100) 


-0.72 


1,2-1,1 


204.63 (100) 


1.38 


2,2-2,2 


277.57 (100) 


-0.82 


0,1-0,0 






0, 1 - 2, 2 






1,0-1, 1 


325.50 (250) 


1.60 


HN'^C: FuFi, F = 0, 1, 1 - 1, 2, 2 at 87090.675 (3) MHz 


2,2,2-1,1,1 


115.5 (50) 


-7.3 


2,2,2-1,2,2 






2,2,1 - 1,1,0 






2,2,1-1,2,1 






2,3,3 - 1,2,2 


158.5 (50) 


5.0 


2,3,2-1, 1, 1 






2,3,2-1,2, 1 






1,1,1-1,1,1 


210.7 (50) 


1.5 



1,2,2-1,2, 1 
1,2,2-1,2,2 



" Numbers in parentheses are one standard deviation in units of the least significant figures. 

* Residuals for DCO* refer to fit with C; kept fixed. The residuals are zero if Cj is released as two pieces of information are 
then used to derive two parameters. 
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Table 4. Derived hyperfine parameters" (kHz) of DCO^ in comparison to previous experimental and present 
quantum-chemicaUy calculated values. 



Parameter 






Value 






This work 


This work 


Caselli & Dore (2005) 


calculated 


eQq 


+ 150.00 (266) 


+ 151.12 (288) 


+ 147.8 (35) 


+ 156.0 


c, 


-0.69* 


-1.12(43) 


-1.59 (78) 


-0.69 



" Numbers in parentheses are one standard deviation in units of the least significant figures. 
* Kept fixed in the analysis. 





if 









7.63055 10 
Rest Frequencji (MHz) 



Fig. 2. Spectrum of the 7 = 1-0 transition of DNC toward the dark cloud LDN 1512, obtained with the IRAM 
30m telescope. Vertical scale is antenna temperature (r^b) in K. Left panel: Gaussian decomposition of the broad 
blended feature into six components which are assumed to have the same width as the single isolated feature. 
Right panel: Spectrum with the isolated feature subtracted and a reduced five-component fit to the remaining broad 
feature. 

Table 5. Hyperfine parameters" (kHz) of DNC in comparison to previous experimental and present quantum- 
chemically calculated values. 



Parameter 




Value 




This work 


Bechtel et al. (2006) 


ab initio 


eQq(N) 


+288.2 (71) 


+294.7(131) 


+309.6 


cm 


+4.91 (63) 


+5.01 (99) 


+5.54 


eQq(D) 


+265.9 (83) 


+261.9(145) 


+257.6 


C,(D) 


-0.60* 




-0.60 


5(ND) 


-1.35* 







" Numbers in parentheses are one standard deviation in units of the least significant figures. 
* Kept fixed to computed values in the analysis. 



Table 6. Hyperfine parameters" (kHz) of HN'^C in comparison to ab initio values. 



Parameter 


Value 






experimental 


ab initio 


eQq(N) 


+272.5 (51) 


+279.5 


C,(N) 


+6.33* 


+6.33 


C,("C ) 


+23.46* 


+23.46 


C;(H) 


-4.49* 


-4.49 


S(HN) 


-8.79* 




5(N"C) 


-1.37* 




5(H'3C) 


-2.98* 





" Numbers in parentheses are one standard deviation in units of the least significant figures. 
* Kept fixed to computed values in the analysis. 
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Table 7. Velocity offsets and intrinsic strengths of hyperflne components of the DCO"^, DNC, and HN"C J = 1-0 
lines. 



Transition 




Component 




Offset 
km s"' 


Relative 
intensity 


DCO* / = 


1-0 


F=0-1 




+0.269 


0.20 






F=2-l 




±0.0 


1.00 






F=l-1 




-0.198 


0.60 


DNC J = 1 


-0 


/, F=2,l-2,2/ 2,1-0,0 




±0.0 


1.00 






/,F=1,1-1,1 




-0.453 


0.99 






/, F=2,3-2,2 




-0.689 


2.33 






/, F= 1,2-1,1 




-0.804 


1.67 






/, F=2,2-2,2 / 0,1-0,0 / 0,1-2,2 




-1.091 


2.67 






/, F= 1,0-1,1 




-1.279 


0.33 


HN"C J = 


1-0 


Fi,F2,F=0,l, 1-1,2,2 




±0.0 


1.00 






Fi, F2,F=2,2,2-1,1,1 / 2,2,2-1,2,2 / 2,2,1- 


-1,1,0/2,2,1-1,1,0 


-0.398 


4.04 






Fi , Fz, F=2,3,3-l,2,2 / 2,3,2-1,1,1 / 2,3,2- 


-1,2,1 


-0.546 


6.63 






Fi,F2,F=l, 1,1-1, 1,1 / 1,2,2-1,2,1 / 1,2,2-1,2,2 


-0.725 


3.66 



